arXiv:1503.07499v2 [cond-mat.mes-hall] 24Jun2015 


Subdiffraction Light Concentration 
by J-Aggregate Nanostructures 


Adriano Cacciola,^ Claudia Triolo,^ Omar Di Stefano,^ Armando Genco,^’^ Marco 
Mazzeo,^’^ Rosalba Saija,"^ Salvatore Patane,"^ and Salvatore Savasta*’*^ 

Dipartimento di Fisica e di Scienze della Terra, Universitd di Messina, Viale F. Stagno 
d^Aleontres 31, F98166 Messina, Italy, CBN, Istituto Italiano di Teenologia, Via Barsanti 
1, 73010 Leeee, Dipartimento di Matematiea e Fisiea ^Fnnio De Giorgi”, UniversitdAL del 
Salento, Via per Arnesano, 73100 Leeee, Italy, and CNR NANOTEC - Istituto di 
Nanoteenologia, Polo di Nanoteenologia e/o Campus Eeotekne, via Monteroni - 73100 

Leeee, Italy 

E-mail: ssavasta@unime.it 

Abstract 

We show, by accurate scattering calculations, that nanostructures obtained from 
thin films of J-aggregate dyes, despite their insulating behavior, are able to concentrate 
the electromagnetic field at optical frequencies like metallic nanoparticles. These results 
promise to widely enlarge the range of plasmonic materials, thus opening new perspec¬ 
tives in nanophotonics. Specifically we investigate ultrathin nanodisks and nanodisk 
dimers which can be obtained by standard nanolitography and nanopatterning tech¬ 
niques. These molecular aggregates display highly attractive nonlinear optical proper- 
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ties which can be exploited for the realization of ultracompact devices for switching light 
by light on the nanoscale without the need of additional nonlinear materials. KEY¬ 
WORDS: J-aggregates, nanophotonics, excitons, surface waves, plasmonics, light har¬ 
vesting 

When light interacts with metal nanoparticles and nanostructures, it can excite collective 
oscillations known as localized surface plasmons (LSPs) which provide the opportunity to 
conhne light to very small dimensions below the diffraction limit This high confinement 
can lead to a striking near-field enhancement which can significantly enhance weak nonlin¬ 
ear processed and enables a great variety of applications such as optical sensing,!^ higher 
efficiency solar cells nanophotonics^^^® including ultracompact lasers and amplifier^^ and 
antennas transmitting and receiving light signals at the nanoscale.®^^ The small mode 
volume of LSP resonances also increases the photonic local density of states (LDOS) close to 
a plasmonic nanoparticle, enabling the modification of the optical properties (decay rate and 
quantum efficiency) of emitters placed in its close proximity (see e.g. Ref.^ and Support¬ 
ing Information Figure SI). The interaction of quantum emitters, as quantum dots or dye 
molecules, with individual metallic nanostructures carries significant potential for the quan¬ 
tum control of light at the nanoscalefirst highlighted by Takahara et al.p^ only 
materials with a negative real part of the dielectric function and moderate losses, are able to 
excite localized surface plasmons and hence to confine light to very small dimensions below 
the diffraction limit. These collective and confined excitations are efficiently supported, de¬ 
spite dissipative losses, by noble metals where the effective response of the electrons can be 
described by a Drude-Lorentz dielectric function whose real part is negative for frequencies 
below the plasma frequency.^^l Also superconductors or graphene have been proposed as a 
platform for surface plasmon polaritons.^^^^ 

Collective oscillations of free electrons are not the only way a negative permittivity 
may arise. It may also occurs in the high energy tail of a strong absorption resonance. 
For example, it has been shown that lattice vibrations in polar dielectric materials can 
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also give origin to negative dielectric permettivity in the far or mid-infrared spectral range 
which can support phonon-polaritons confined to the surfaceIt has also been shown 
that in the near-infrared and optical spectral ranges, sharp excitonic resonances in organic 
molecules can give origin to spectral regions where the dielectric permettivity is negative 
The response of Lorentzian/excitonic nanoshells with optically inactive core and embedding 
medium have been also studied in the quasistatic limit and compared with the proper¬ 
ties of the Drude/plasmonics nanoshells.®! Recently, it has been theoretically shown that 
a nanosphere of TDBC-doped PVA displays enhanced optical fields and subwavelength 
field confinement analogously to plasmonic nanoparticles.®! ^t present, however, spheri¬ 
cal nanoparticles of J-aggregates have not been realized and the synthesis of high quality 
nanospheres does not seem readily achievable. 

Here we propose more realistic J-aggregates nanostructures suitable for obtaining on chip 
strong subwavelength light confinement. By accurate scattering calculations, we show that 
thin films of J-aggregate dyes®! can be exploited for the realization of plasmon-like nanos¬ 
tructures working at optical frequencies. Such nanostructures can be realized after the film 
growth by standard nanolitography and nanopatterning techniques. The proposed struc¬ 
tures consist of ultrathin nanodisks on a silica substrate. Moreover, we show that these 
molecular aggregates can also be exploited for the realization of more complex plasmon-like 
nanostructures to achive very high field concentration. Specifically we find that nanodisk 
dimers with radius a = 30 nm are able to generate near-field enhancements \E/Eq\^ beyond 
350. The here proposed J-aggregate nanostructures are promising towards the emerging of a 
low-cost plasmonic technology for the development of ultracompact photonic and optoelec¬ 
tronic devices. Moreover, J-aggregates display significant optical nonlinearitie^®®! which 
can induce switchable localized surface resonances (LSRs) required for the realization of 
plasmonic devices!^ without the need to exploit the coupling with additional nonlinear opti¬ 
cal materials.®®! This target has not been achieved by metallic nanoparticles owing to their 
intrinsically weak optical nonlinearities which can be overcome only by coupling LSPs to 
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nonlinear optical resonances 

J-aggregates self-assembled nanostructures have been realized in various forms including 
quasi-one dimensional chains, single monolayers, tubes, etc. which could be employed as 
molecular plasmonic materials.® One limitation of these surface resonances is that they 
can manifest in a limited spectral region in the high energy side of the excitonic resonance. 
However, because of the huge variety of organic dyes that can aggregate, it could be possible 
to create J-aggregates LSRs within an arbitrary part of visible and near IR spectrum® and 
there is room for chemically engineering and optimizing their plasmonic behavior in order to 
minimize losses. Moreover J-aggregates emit light, generally with good quantum efficiencies, 
thus their LSRs could be exploited for the realization of efficient light-emitting nanoantennas. 

These J-aggregates possess remarkable morphological and optical properties that make 
them suited for use in photonic and optoelectronic applications.®® The broad absorption 
line corresponding to the single molecule excitation is shifted to the red side of the spec¬ 
trum and transforms into a collective narrow-linewidth optical transition possessing a giant 
oscillator strength derived from the coherent coupling of the aggregated molecules.® Layer 
by layer growth of cyanine-dyes (TDBC) aggregates to form an ordered thin film has been 
demonstrated with a peak absorption constant of a ~ 1 x 10 ® cm“^, one of the highest 
reported for a thin film.® The Frenkel excitons of these high absorptive materials have been 
exploited in demonstrations of strong coupling between excitons and microcavity photons,® 
room-temperature operation of exciton-optoelectronic devices,® strong coupling between 
excitons and surface plasmon-polaritons® or LSP resonances.® 

The dielectric function at frequencies around sharp excitonic resonances is usually quite 
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well reproduced by the Lorentz model: e{u)) = £oo + where ujq is the resonance 

frequency, 70 the width of the resonance, / the dimensionless oscillator strength and £00 fhe 
high-frequency background constant. In contrast to the Drude dielectric function, whose 
real part is negative for frequencies below the plasma frequency, the real part of the Lorentz 
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Figure 1: Plasmon-like properties in the blue-side spectral region of strong excitonic reso¬ 
nances. (a) Real and imaginary part of the dielectric function directly derived by the complex 
refractive index (h = n + in) obtained via a Kramers-Kronig regression of reflectance data^^ 
of layer by layer assembled TDBC/PDAC thin films. The real part of the dielectric function 
is negative in the blue-side region, a key feature for confining light below the diffraction 
limit. The inset shows the molecular structure of the polycation PDAC and of the anyon 
TDBC, a J-aggregate forming cyanine dye. (b) Colormap of the peak extinction efficiency 
Qext = Calculated in the quasistatic approximation for a 40 nm sphere in vacuum 

by using the Lorentz dielectric function as a function of the normalized oscillator strength 
/ and the high-frequency dielectric constant £oo with fixed Qo = 80. The white line in 
the contourplot separates the region where £min < — 2 (on the right side) and plasmon-like 
effects starts to be expected, from the left region where £min > — 2. This model map points 
out the plasmonic potentiality of excitonic resonances with high oscillator strengths. The 
corresponding values of Qext calculated for Ag, Au, and TDBC/PDAC nanoparticles are 
indicated in the color-bar for reference. 
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dielectric function is generally positive for almost all excitonic absorption lines. It takes its 
minimum ~ £oo — /<^o/(27o) at Wmin — <^o + 7o/2. However there is no specific constraint 
preventing negative values. For narrow absorption lines and high oscillator strengths, the 
real part of the Lorentz dielectric function can acquire negative values on the high energy 
side of the excitonic resonance if /(5o/(2£oo) > 1, where Qo = <^o/7o is the quality factor of 
the excitonic line. The giant oscillator strength and the strong spectral narrowing of TDBC 
aggregates, both induced by the coherent intermolecular coupling, result in ordered thin 
films with a negative dielectric function in the high energy tail of the excitonic resonance 
with an astonishingly low minimum value £min — —14. Figure la displays the dielectric 
function directly derived by the the real and imaginary components of the refractive index 
(h = n + in) obtained via a Kramers-Kronig regression of neat films reflectance data.^^ 
The obtained optical constants accurately reproduce both the reflectance and transmittance 
data for a 5 nm thick thin film and the reflectance of more complex planar optical structures 
containing a thin TDBC/PDAC film.l^ 

The obtained TDBC/PDAC dielectric function shown in Fig. la qualitatively resembles a 
Lorentz dielectric function. However its imaginary part displays a tail in the blue part of the 
spectrum due to the presence of the vibrational degrees of freedom and the Kramers-Kronig 
related real part is also asymmetric. 

Results 

In order to explore the potentialities of J-aggregates as plasmonic materials, we start with 
the most convenient geometry for an analytical treatment: a spherical nanoparticle of radius 
a much smaller than the effective wavelength X/^/ed, where £d is the dielectric constant 
(assumed to be real) of the surrounding medium. Scattering calculations on more feasible J- 
aggregates nanostructures will be presented later. The signature of LSP effects can be easily 
recognized by the polarizability cCp = P/(£o£dF'o), defined as the ratio between the dipole 
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moment P induced in the nanoparticle by an incoming wave and the amplitude of the incident 
displacement field Eq. In the quasi-static approximation, for a spherical nanoparticle, Op = 
47ro^(£(a;) — £d)/{s{u) -|- 2sd)- The nanoparticle polarizability displays a dipole resonance 
when £'{oi>) = —2£d known as Frolich condition. The optical theorem provides a direct link 
between the polarizability and the extinction cross section (Text = \/^^o Ini{a;p}, where 
ko = 2 tt/\. This resonance also produces a resonant enhancement of the dipolar field and of 
the optical LDOS around the nanoparticle. Many of the attractive plasmonic applications 
of metal nanoparticles and nanostructures rely on this field-enhancement at the plasmon 
resonance por a TDBC nanoparticle in vacuum, the dipole LSR occurs at A ~ 570 
nm where e = —2 -|- i3.5. In Fig. lb a model-map points out the plasmonic potentiality 
of strong excitonic resonances. The figure displays the peak extinction efficiency Qext = 
crext/(7ra^) for a 40 nm sphere calculated within the quasistatic approximation by using the 
Lorentz dielectric function versus the normalized oscillator strength / and the high-frequency 
dielectric constant e^o with fixed Qo = 80. 

Figure 2a displays the extinction efficiency Qext for TDBC spherical nanoparticles of 
different radii a, calculated beyond the quasistatic approximation, by employing the Mie 
theory implemented within the T-matrix formalism^El (see also the Method Sect.). It is 
remarkable that a Qext slightly larger than 3 is observed for a nanoparticle with a radius 
of 60 nm. This means that the interception area that the spherical particle presents to the 
incident radiation is about three times larger than its actual size. Such extinction efficiency 
is lower than that for a gold nanoparticle with the same radius (see also the colorbar in 
Fig. lb), but it is significantly larger than the efficiency of dielectric particles with positive 
dielectric functions. It is worth noticing that the extinction spectra displays peaks at the 
wavelength corresponding to the LSR which is significantly blue-shifted with respect to the 
excitonic absorption maximum (compare Fig. 1 and Fig. 2a). Figure 2a also displays a 
significantly smaller peak corresponding to the excitonic resonance. Figure 2b shows the 
calculated near-field enhancement spectra \E/Eq\^ for the TDBC nanoparticles. The field 
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Figure 2: Scattering calculations on TDBC/PDAC nanosptieres. (a) Extinction efficiency 
Qext for TDBC spherical nanoparticles in vacuum of different radii a, calculated beyond the 
quasistatic approximation, by employing the Mie theory implemented within the T-matrix 
formalism, (b) Near-held enhancement spectra jE/Eop for the TDBC spherical nanoparti¬ 
cles. (c) Near-held enhancement distribution on the plane containing the polarization axis 
and normal to the propagation direction for a TDBC nanoparticle with a radius a = 40 nm 
illuminated by a plane-wave at A = 578 nm. (d) Near-held enhancement spectra as in (b) 
but with different dielectric constants of the external surrounding medium. 













has been calculated at a distance d = 1 nm from the surface of the nanoparticle at an 
angle between the incident electric held and the position vector (with the origin at the 
nanoparticle center) 6 = Q. A signihcant plasmon-like near-held enhancement (beyond one 
order of magnitude for the a = 60 nm nanoparticle) can be observed. Such enhancement 
is about one half of the enhancement induced by a gold nanoparticle of the same size. 
The wavelength where the maximum enhancement occurs is red-shifted with respect to the 
corresponding extinction-spectrum peak, which is a typical plasmonic feature.^Sl Figure 2c 
shows the near-held enhancement distribution on the plane containing the polarization axis 
and normal to the propagation direction for a TDBC nanoparticle with a radius a = 40 nm. 
Calculations have been performed considering an incident plane-wave at A = 578 nm. The 
inhuence of the dielectric environment on the LSR resonance is shown in Fig. 2d. A shift of 
the resonance due to the modihcation of the Frolich condition can be observed. By conhning 
light using LSPs, it is possible to signihcantly alter the photonic LDOS, thus enabling a 
remarkable enhancement of the spontaneous emission rate of emitters placed in the near- 
held of a metallic nanoparticle. For an emitter oriented towards a 40 nm TDBC nanoparticle 
center, located at a distance d = 5 nm with a transition wavelength Aq = 578 nm, we hnd a 
spontaneous emission enhancement F/Fq — 10 (Supporting Information Figure SIl). 

When two or more plasmon-resonant nanoparticles are closely spaced, LSPs of individual 
particles interact, giving rise to a strong optical held enhancement in the gap regions between 
the particles.®® As a further signature of the plasmonic behaviour of TDBC nanoparticles, 
we show that such a strong enhancement can be observed with J-aggregate TDBC dimers. 
The near-held enhancement spectrum \E/Eq\^ calculated at the center of the gap between two 
nanospheres with a = 40 nm is reported in Fig. 3a. All the presented calculations for dimers 
have been obtained considering a longitudinally polarized incident held (along the dimer 
axis). The peak held enhancement at the center of the gap between two TDBC spherical 
nanoparticles (a = 40 nm) as a function of the gap is reported in Fig. 3b. Figure 3c shows 
the peak held-enhancement at the center of a 5 nm gap as a function of the nanoparticles 
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Figure 3: Scattering calculations for a dimer of two nanospheres in vacuum illuminated by 
longitudinally polarized plane waves, (a) Near-held enhancement spectrum calculated at 
the center of a 5 nm gap between two nanospheres (radius a = 40 nm). (b) Peak held 
enhancement at the center of the gap as a function of the gap (o = 40 nm). (c) Peak 
held-enhancement at the center of a 5 nm gap as a function of the nanoparticles radius, (d) 
Near-held enhancement distribution around the nanostructure on the plane z = h/2 for h=5 
nm calculated for an incident plane wave at A= 581 nm. 
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radius. The near-field enhancement distribution on the plane containing the dimer axis 
and orthogonal to the propagation direction is displayed in Fig. 3d. At present, spherical 
nanoparticles of J-aggregates have not been realized and we are not aware of any present 
growing technique. Moreover, the dielectric function here considered has been derived by 
means of optical measurements on an ordered TDBC thin film grown layer by layer on a 
substrate. 


We thus study more feasible J-aggregates nanostructures, showing a plasmon-like behav¬ 
ior analogous to the above described nano-spheres. Such nanostructures can be realized after 
the film growth by standard nanolitography and nanopatterning techniques.®' The proposed 


structures consist of ultrathin nanodisks and nanodisk dimers grown on a silica substrate. 
Scattering calculations have been performed by the finite element method (see Methods). 
Figure 4a displays the extinction efficiency Qext spectra for three single nanodisks of radius 



Figure 4: Scattering calculations on single nanodisks, (a) Extinction efficiency Qext spectra 
for three single nanodisks of radius a = 30 nm and heights h = 5, 10, 15 nm, calculated under 
normal illumination, (b) Near-field enhancement calculated at .2 = h/2 and at a distance 
d = Inm from the disk lateral surface along the polarization direction (0 = 0). (c) Near-field 
enhancement distribution on the plane z = /i/2 and normal to the propagation direction for 
the TDBC nanodisk with h = 5 nm, under plane wave illumination at A = 591 nm. 


a = 30 nm and heights h = 5, 10, 15 nm, calculated under normal illumination. Increasing 
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the thickness h, a blue-shift of the LSR is observed. This blue-shift at increasing ratios 
h/a is a typical plasmonic feature reported for metal nanodisks The small shoulder on 
the red side of the peak is a signature of the exciton absorption peak at about A = 595 
nm corresponding to the wavelength where the imaginary part of the TDBC/PDAC dielec¬ 
tric function is maximum (see Fig. la). The near-field enhancement \E/Eq\^ calculated at 
an angle between the incident electric field and the position vector (with the origin at the 
nanodisks center) 9 = 0, at z = h/2, and at a distance d = Inm from the disk lateral 
surface, is shown in Fig. 4b. Figure 4c displays the near-field enhancement distribution on 
the plane containing the polarization axis and normal to the propagation direction for the 
TDBC nanodisk with h = 5 nm. Calculations have been performed considering an incident 
plane-wave at the near-field enhancement peak in Fig. 4b. A strong field enhancement be¬ 
yond 25 is observed. These results clearly show that realistic J-aggregate nanostructures 
are able to support LSP-like resonances. For comparison, scattering calculations for gold 
nanodisks have been included in the Supporting Information (Fig. SI2). 

It is worth observing that the dielectric function used for all the presented calculations 
has been obtained from optical data fitted by an isotropic model, thus any optical anisotropy 
effect in neglected. With J-aggregates having an intrinsically anisotropic structure ,^^211] 
anisotropy of the dielectric response can be expected in thin films as well. At our knoledge 
there are not specific data on TDBC films, however recently, by using spectroscopic ellip- 
sometry and polarized IR spectroscopy, a strongly anisotropic behavior of thiacarbocyanine 
dye aggregates films was revealed.^ It was quantified by means of two dielectric functions: 
one for the xy plane parallel to the film and the other for the 2 : direction perpendicular to it. 
These data do not show any optical anisotropy for polarizations along the xy plane. In the 
Supporting Information (Figure SIS) we show optical absorbance measurements at normal 
incidence for two different orthogonal linear polarizations. The data confirm the absence of 
optical anisotropy on the xy plane. All the calculations on nanodisks have been performed 
by considering input light at normal incidence and hence polarized along the xy plane. In 
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Figure 5: Scattering calculations under plane wave illumination with normal incidence and 
longitudinal polarization (along the dimer axis) on TDBC/PDAC nanodisk dimers on a 
glass substrate. The two disks have a radius a = 30 nm and are separated by a gap of 4 nm. 
(a) Extinction efficiency for dimers of two different heights h. (b) Near-held enhancement 
spectra calculated at the gap center and at z = h/2. (c) Near-held enhancement distribution 
around the nanostructure on the plane z = h/2 for h = 5 nm and A = 591 nm. (d) Near-held 
distribution calculated on the substrate surface z = 0 for h = 5 nm and A = 591 nm. 
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this case the near field generated by the nanodisk along the z axis is very small (beyond 
one order of magnitude) and thus the presence of a different dielectric function along the 
z axis does not affect significantly the obtained results. The lack of optical anisotropy in 
the xy thin film response may be due to the local disorder during the growth process, as 
evidenced in AFM images in Ref.^ We also observe that, from available AFM images (see 
e.g. Ref.l^, it is possible to infer some degree of local order in the 100 nm range. Hence it 
is possible that TDBC nanodisks (or more generally nanostructures) display some degree of 
optical anisotropy on the xy plane, depending on the specific arrangement of the monomers. 
In this case the dielectric function that we used for all the calculations can be interpreted 
as an average quantity. For example for a brickwork type arrangement of the monomers 
in the aggregates, a larger dipole moment can be expected along the head-tail direction of 
the aggregates. For the nanodisks here investigated, we would expect an increase of the 
near-field enhancement for input light polarized along the head-tail direction and a corre¬ 
sponding decrease along the orthogonal direction. Another interesting issue is the possible 
dependence of the dielectric function on the size of the nanostructures. As the sharp peak 
in the dielectric function and the corresponding negative permettivity arises from the in¬ 
teraction of the dye molecules, the dielectric function might strongly depend on the size of 
the nanostructures. According to accurate theoretical calculations including both the elec¬ 
tronic and vibrational degrees of freedom, a chain with 15 monomers fully reproduces the 
observed shifted and sharp absorption peaks of J-aggregates.^^ These theoretical analysis 
is confirmed by femtosecond nonlinear optical experiments on TDBC J-aggregates, showing 
that the exciton delocalization length at room temperature corresponds to 16 molecules.l^ 
As a consequence, nanostrutures with linear dimensions beyond 30 nm should display the 
desired optical properties. 

Scattering calculations on nanodisk dimers are shown in Fig. 5. The two disks have a 
radius a = 30 nm and are separated by a gap of 4 nm. We consider dimers of two different 
heights h = 5, 10 nm. The extinction efficiency and near-field enhancement spectra are shown 
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in Figs. 5a and b under plane wave illumination with normal incidence and longitudinal 
polarization (along the dimer axis). The enhancement spectra have been calculated at the 
gap center and z = h/2. The near-held enhancement distribution around the nanostructure 
on the plane z = h/2 under the above described conditions is displayed in Fig. 5c. Figure 5d 
shows the resulting held distribution calculated on the substrate surface (z = 0). A further 
concentration of the held around the opposite disks surface can be observed. Analogous 
calculations based on an excitonic resonance providing a dielectric function with positive real 
part^^ displays a negligible near-held enhancement (Supporting Information Figure SI4). 


Conclusions 

The presented results demonstrate that J-aggregates nanostructures display LSRs able to 
conhne light below the diffraction limit, resulting into a strong near held enhancement com¬ 
parable to that induced by gold nanoparticles. These results, in view also of the huge variety 
of organic dyes that can aggregate and their chemical hexibility, hold the promise to greatly 
enlarge the availability of new plasmonic materials with different properties with respect 
to noble metals. For example J-aggregates display attractive nonlinear optical properties 
which could be exploited for the realization of switchable LSRs for nanophotonic devices. 
Specihcally, it has been shown that an ultrafast optical pump is able to signihcantly affect 
the dielectric function of J-aggregates,*^*^ ^^ * ^^ * hence we expect light-controlled modihcations 
of the LSRs, very sensitive to the real part of the dielectric function, which can be exploited 
for the realization of ultracompact all-optical switches. Since analogous molecular aggre¬ 
gates can be found in nature,^ acting as antennas in photosynthetic complexes, the results 
here presented open the intriguing question if photosynthetic organisms!^ exploit LSRs to 
improve light harvesting efhciency.* ^^ * ^^ * At the same time J-aggregate nanoantennas could 
represent a low-cost solution for improving photovoltaic devices and for the realization of 
efficient artificial photosynthetic light-harvesting antennas.® 
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Methods 


Scattering calculations on spherical nanoparticles and dimers have been performed within 
the transition-matrix method, initially introduced by Waterman!^ as a general technique 
for computing the e.m. scattering based on the Huygens principle. In the last decades this 
method has been significantly improved by expanding the incident, the internal and the 
scattered electromagnetic fields in vector spherical harmonics (VHS).® Thanks to the lin¬ 
earity of MaxwellaAZs equations and boundary conditions, when the scatterers, either single 
or aggregate, are constituted by spherical monomers, the expansion of the electromagnetic 
fields in terms of VSH allows to analytically relate the multipolar amplitude of incident field 
to those of the scattered field by means of the T-matrix. The Transition matrix contains 
all the information on the microphysical properties of the scatterer, being independent from 
the state of polarization of incidence field, and from the incident and observation direction. 
The elements of T-matrix define analytically in far-field, the optical cross section, and in the 
near-field zone the intensity distribution of the total and scattered field. From the computa¬ 
tional point of view, the numerical calculation of the transition matrix requires the inversion 
of a matrix whose order is, in principle, infinite. However, it is possible to truncate the 
multipole expansion of the fields, that determine the size of the array, in order to ensure the 
numerical stability of the results. 

Finite element method (FEM) is a differential equation method that computes the scat¬ 
tered time-harmonic electromagnetic field by solving numerically the vector Helmholtz equa¬ 
tion subject to boundary conditions at the particle surface. Thanks to the discretization 
process (mesh construction) of the domains in smaller parts, this technique is suitable to 
describe systems with complex geometry. The structure, placed on a glass substrate is 
embedded in a 3D-finite computational domain that is discretized into many small-volume 
tetrahedral cells. In the far-field zone, at the outer boundary of the finite computational 
domain, approximate absorbing boundary conditions are imposed. In this way it is possible 
both to suppress wave reflections back into the domain and to allow the propagation of the 
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numerical outgoing waves as if the domain were infinite.!^ In the near-field zone, for the 
3D tetrahedral elements, we choose an extremely fine size (up to 0.1 nm) to better describe 
the rapid changes in the electromagnetic field. By solving 3D Helmholtz’s equations to¬ 
gether with the boundary conditions in each element of mesh, we can reconstruct the optical 
behavior of the system, the near-field enhancement distribution around the nanostructures 
(see figures 4b, 5c-d) and the extinction spectra in far-field. Scattering calculations were 
validated by comparing the results with the simulation of the same spherical nanoparticles 
and dimers studied above with the Mie theory. 
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